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ABSTRACT 

This paper presents new observations of the planet-hosting, visual binary GJ 86 
(HR637) using the Hubble Space Telescope. Ultraviolet and optical imaging with 
WFC3 confirms the stellar companion is a degenerate star and indicates the binary 
semimajor axis is larger than previous estimates, with a > 28 AU. Optical STIS spec- 
troscopy of the secondary reveals a helium-rich white dwarf with C2 absorption bands 
and T c ff = 8180 K, thus making the binary system rather similar to Procyon. Based on 
the 10.8 pc distance, the companion has 0.59 Mq and descended from a main-sequence 
A star of 1.9 Mq with an original orbital separation a > 14 AU. If the giant planet is 
coplanar with the binary, the mass of GJ86Ab is between 4.4 and 4.7 Mj up . 

The similarity of GJ 86 and Procyon prompted a re-analysis of the white dwarf 
in the latter system, with the tentative conclusion that Procyon hosts a planetesimal 
population. The periastron distance in Procyon is 20% smaller than in aCenAB, 
but the metal-enriched atmosphere of Procyon B indicates that the planet formation 
process minimally attained 25 km bodies, if not small planets as in a Cen. 

Key words: binaries: visual — stars: individual (GJ86A, GJ86B) — planetary 
systems — white dwarfs 



, 1 INTRODUCTION 

About 20% of known extrasolar planets orbit one compo- 
nent of a stellar binary l|Raghavan et al . 2006; Haghig hipour! 
2006). Both observation and theory indicate that the major- 
ity of these planets in binaries are similar to those orbiting 
single stars, owing to the wide (a > 100 AU) stellar separa- 
tions (|Desidera fc BarbierilfeoOTtlHaghighipourlbOOel h How- 
ever, there are a growing number of systems which present 
a major challenge to planet-formation modelers, because 
their giant planets orbit one component within closer bi- 
nary systems. In order of discovery, these planetary sys- 
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terns are GJ86 (HR637 HP 13445) (|Queloz et al] l200Ch 
7 Cep hei (lHatzes et all l2003h. HP 410 04 (|Zucker et all 
120041) . HP 196885 jCorreia et all l2008h . and HP 17605 
1 Muterspaugh et al.l |2010| ). All five of these binaries are 
thought to have a < 30 AU and a giant planet orbiting 
the primary star. However, only in GJ 86 is the stellar sec- 
ondary a white dwarf rather than a second main-sequence 
star, and the Jovian planet thus orbits the originally less 
massive component. Furthermore, the initial binary separa- 
tion was smaller, making GJ 86 a challenging environment 
in which to form planets. 

Our view of planet formation in binaries has changed 
significantly in the past 17 years. Simulations of binary 
stars with a < 40 AU failed to retain sufficient circum- 
stellar material to form planets around one star, and thus 
precluded planet birth by either the core accretion or disk- 
instability mechanisms ^Nelson 200TJ; lArtvmowicz fc Lubowl 
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Table 1. HST Cycle 19 Observations of GJ86B 



Table 2. Flux Densities & Effective Wavelengths for GJ 86B a 
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1 1994 see iPrato fc Weinberger! l20ld for a detailed review). 
Yet observers have imaged favorable planet-forming envi- 
ronments around the components of close binary stars, im- 
plying planet formation in t hese systems may b e as com- 
mon as around single stars |Akeson et al l Il998l ; iMathieul 
11994 ). For example, the two well-separated disks in the bi- 
nary system L1551 (Rod riguez et alj |l998) indicate that, 
despite disk truncation, it is still possible for both compo- 
nents to retain a significant amount of their original cir- 
cumstellar material (0.03 — 0.06 A/q) in disks with consid- 
erable radii (r ~ 10 AU). These disk masses are comparable 
to the minimum-mass model of the p rimordial Solar nebula 
|Havashilll98ll ; IWeidenschillinglll977l ) suggesting that planet 
formation in these environments may proceed as in disks 
around single stars. 

The Jovian planelQ orbiting GJ 86A (K0 V) was discov- 
ered via precision radial- velocity monitoring |Queloz et al.l 
2000), revealing a 15.8 d orbit (0.11 AU separation) with 
msini = 4Mj up . At the time of discovery, an additional 
long-term radial velocity drift was reported, suggesting the 
presence of a more distant, unseen stellar companion. Sub- 
sequently, the distant companion was directly detected at 
a projected separation of l'/7, but 9 ma g fainter than the 
primary in the K band (lEls et all 1200 ll). The companion , 
GJ 86B, was later imaged by lMugrauer fc Neuhauserl (2005) 
and lLagrange et al.l (|2006l ) whose data clearly reveal orbital 
motion; between 2000 and 2005 the angular separation of 
the pair changed by 0'.'25 and the position angle changed 
by 17°. GJ 86B cannot be a low-mass star or brown dwarf: 
a substellar mass cannot be reconciled with the observed 
radial-velocity drift of GJ 86A, and follow-up methane-band 
photometry showed GJ 86B to have a near-IR color index 
near zero. Moreover, A"-band spectroscopy of GJ 86B re- 
veals an essentially Rayleigh- Jeans co ntinuum source. Thus , 
GJ 86B can only be a white dwarf l|Lagrange et alj 120061 ; 
iMugrauer fc Neuhauseill2005l ). 

In order to ideally constrain the formation and evolu- 
tion of the giant planet orbiting GJ 86A, it is necessary to 
obtain the orbit and component masses of the binary sys- 
tem. The 3300:1 contrast ratio between components in the 

1 The planet is often referred to incorrectly as GJ86b, while the 
correct designation is GJ 86Ab 
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a The flux of GJ86A was measured to be 31.4 mjy in F225W. 
b Effective wavelengths were derived by convolving the adopted 
stellar model in F\ with the WFC3 filter transmission curves. 



infrared makes this difficult but possible with adaptive op- 
tics, yet white dwarfs do not exhibit spectral features in this 
wavelength region. Thus an optical study of GJ 86B is neces- 
sary and a Hubble Space Telescope (HST) program was ini- 
tiated to achieve this goal. In §2 are described the imaging 
and spectroscopic observations, while §3 describes the at- 
mospheric modeling, derivation of current (and progenitor) 
stellar parameters and binary orbit constraints. New static 
limits on the GJ 86 planetary system are calculated based 
on the HST observations and possible similarities with the 
Procyon system are discussed in §4. 



2 OBSERVATIONS & DATA REDUCTION 

Due to orbital motion observed between 2000 and 2005 
l|Lagrange et alj|200a : IMugrauer fc Neuhau scr 2005), it was 
expected that precise spatial offsets might be necessary for 
proper placement of the faint companion within the spec- 
troscopic aperture. Thus, WFC3 imaging was executed for 
both photometric and astrometric purposes, prior to STIS 
spectroscopy. 

2.1 WFC3 Imaging 

GJ 86 was imaged on 2012 March 31 using the Wide Field 
Camera 3 (WFC3) in the UVIS channel with a 512 x 512 
(UVIS2-C512C-SUB) subarray. The imaging sequence took 
place over a single orbit and began with a series of 1 s ex- 
posures in the F225W filter. In this short exposure time, 
the flux from the V = 6.2 mag, T c s ~ 5200 K primary was 
expected to leave an unsaturated, linear response on the 
CCD. This sequence employed the default UVIS dither pat- 
tern with four points in a parallelogram, and was repeated 
four times for a total of 16 exposures. Both primary and sec- 
ondary star were well-detected in this first set of exposures, 
and in the linear response regime with high signal-to-noise 
(S/N) > 500. This set of images was used to robustly mea- 
sure the separation of the binary, where the secondary star 
was found at offset: 

a = 2'.'351(2) 6 = 88?96(4) (J2000) Date: 2012.2468 

Immediately following the short exposure sequence, a 
single dither pattern of four deeper exposures was executed 
in each of eight filters (see Table [T] and Figure [TJ. At the 
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Figure 1. Images of GJ 86 in the eight specified ultraviolet and optical filters using WFC3 UVIS. The images are oriented so that 
upward is 17? 60 east of north, and each frame is 8" on a side. The bottom-right image is a composite RGB frame using the F225W 
(blue), F275W (green), and F336W (red) filters. A series of short exposures at F225W were obtained with an unsaturated primary, and 
these data were used for the relative astrometry of the pair. 



phase II design stage, the brightness of the companion was 
poorly constrained and the longer exposures were to insure 
high S/N photometry (up to the limit imposed by the PSF 
wings of the primary) at all wavelengths. Fortunately, even 
without stray light correction, the dither-combined image 
sets yielded S/N > 500 for GJ 86B shortward of 5500 A, and 
S/N > 200 beyond. Thus, the relatively faint companion was 
mostly unaffected by the primary. 

Photometry was performed on multi-drizzled images 
produced by the STScI software and data pipeline OPUS 
2012.1 and CALWF3 2.6.2. The flux of GJ 86B was mea- 



sured in an r = 7 pixel aperture radius and corrected using 
the encircled energy dat a available in the W FC3 Instrument 
Handbook version 3.0 (|Dressel et alJl201Ch . The aperture 
corrections were derived by interpolating the tabular values 
at the effective wavelength of the best fitting atmospheric 
model (see §3) in each bandpass, and the process was iter- 
ated until there were no significant changes in the photom- 
etry or stellar model parameters; in practice this occurred 
after a single repetition. The background flux of GJ 86A was 
subtracted by examining points opposite, in x and y on the 
subarray, from the bright primary. This ambient light cor- 
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Figure 2. Combined-frame STIS CCD spectral images of GJ 86B taken in three gratings, with wavelength increasing to the right. The 
ambient light of the primary is relatively mild, even in red light. Despite rejection routines that drastically reduced the number of cosmic 
rays, hot and cold pixels present in the single exposures, some artifacts remain in the combined 2D images, especially at the shortest 
wavelength setting. 



rection resulted in flux measurement changes less than 1.2% 
for the five filters below 5000 A, but rose to around 4% for 
the three longer wavelength filters. The measured fluxes for 
GJ86B are listed in Tabled 



2.2 STIS Spectroscopy 

Observations with the Space Telescope Imaging Spectro- 
graph (STIS) were executed on 2012 May 29, beginning with 
an acquisition sequence of GJ 86A, followed by an offset and 
acquisition/peak-up of GJ86B based on the WFC3 astrom- 
etry. Examination of the acquisition and peak-up images 
shows they were successful. Spectroscopy was performed 
with the G430L, G430M, and G750M gratings as detailed 
in Table [TJ using default STIS dither patterns along the slit, 
and achieving average resolutions of 6.8, 0.7, and 1.4 A, re- 
spectively. 

The individual spectral images were first shifted onto a 
common frame as indicated by the image header astrome- 
try and using integer pixel shifts. Next, the shifted images 
were combined into a single 2D frame using mscombine in 
the IRAF STSDAS software package. Prior to this step, the 
spectral images contained a significant number of bad pixels 
and cosmic rays, especially in the G430M observations. Re- 
jection of these unwanted pixels was performed using CCD- 
CLIP at 4a, although this could not completely remove the 
accumulated artifacts. 

Figure [2] shows the final, combined spectral images for 
each of the grating settings. Although stray light from the 
primary is apparent in all the images, and most readily seen 
in the red, the signal from GJ 86B was essentially unaffected. 
The contrast between the peak signals from the science tar- 
get and the primary diffraction spikes was 17:1 at the lowest 
in the G750M observations, and at least twice as high in the 
region of the science spectrum. 

Lastly, ID spectra were extracted using xlD, and ap- 
propriate sky subtraction to correct for the relatively mild 



background. The flux-calibrated and normalized STIS spec- 
tra are displayed in Figures [3] and [4] The S/N per pixel was 
determined from continuum regions within each of the spec- 
tra to be close to 85, 88, and 53 in the G430L, G430M, and 
G750M gratings respectively. 



3 STELLAR & BINARY PARAMETERS 

The upper spectrum shown in Figure 0] displays the distinc- 
tive and broad C2 Swan bands, indicating that GJ 86B has 
a DQ spectral classification. This class of white dwarfs is 
characterized by helium-dominated at mospheres with trac e 
carbon, likely dredged from the core (|Pelletier et alJll986h . 
Neither Can K nor Ha were detected in the relatively deep 
and higher resolution G430M and G750M spectra. Follow- 
ing convention the star has the wh ite dwarf designation 
WD 0208-510 (|McCook fc Sionlll999l '). 

3.1 Atmospheric Modeling 

Stellar parameters for GJ 86B were derived in two ways, us- 
ing model atmosphere g rids appropriate for DQ white dwarfs 
fsee lDufour et al"]|2005l for details on the models and fitting 
techniques). First, the effective temperature and solid an- 
gle were estimated from fitting the broad baseline WFC3 
photometry. These parameters were t hen combined with 
the Hipparcos parall ax of 92.74 mas l|van Leeuwenl 120071 ; 
iPerrvman et al.lll997l ') to determine the stellar radius, which 
was then con verted into mass usin g theoretical mass-radius 
relationships (|Fontaine et al.1 120011 ). However, the shape of 
the energy distribution - and thus the derived effective tem- 
perature - is slightly sensitive to the amount of photospheric 
carbon. Model fits to the C2 Swan bands in the G430L 
STIS spectrum were used to measure the carbon abundance, 
with T e fT and log g fixed at the values obtained from the 
WFC3 photometry. The effective temperature and solid an- 
gle were re-derived using C/He fixed to the spectroscopic 
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Figure 3. Flux-calibrated WFC3 photometry and STIS spectra 
of GJ 86B for each of the configurations specified in Table[T] Fig- 
ure [4] plots the spectra in greater detail. 



Table 3. Current & Progenitor Stellar Parameters for GJ86AB 



GJ 86B / WD 0208-510 



White Dwarf Remnant: 
SpT 

V cf f (mag) 

T efi (K) 

log g (cms -2 ) 

Mass (Mq) 

Radius (Rq) 

Cooling Age (Gyr) 

log(C/He) 

log(H/He) 

log (Ca/He) 

Main-Sequence Progenitor: 
SpT 

Mass (Mq) 
Lifetime (Gyr) 



DQ6 
13.2 
8180 ± 120 
8.02 ±0.02 
0.59 ±0.01 
0.01245 ±0.0015 
1.25 ±0.05 
-4.8 ±0.2 
< -4.3 
< -11.8 



A5V 
1.9 ±0.1 
1.4 ±0.2 



GJ 86A / HR637 / HD 13445 



SpT 

V (mag) 
T eff (K) 

Mass (Mq) 
Age (Gyr) 



K0V 
6.1 
5200 
0.80 
2.5 



Note. Parameters listed 
tive of those found 



for GJ86A are 
the literature 



values repr e senta- 



llGhezzi et al.1 

van Belle & von Braunl |2009|; iMamaiek fc Hillenbrand! 
Valenti fc Fischer! 120051: ISaffe et al.1 |2005| ; iRibas etldJ, 
Flvnn fc Morellll 19971 ). 



2010 



20O 



2003 



value. This procedure was iterated until all parameters con- 
verged, with a good fit obtained at T e s = 8420 ± 90 K, 
log (C/He) -4.58 ±0.15. 

Second, because the STIS data are flux-calibrated with 
a high degree of confidence, the effective temperature and 
carbon abundance can be simultaneously obtained from 
the G430L spectrum alone, yielding T cff = 8145 ± 180 K, 
log (C/He) = — 4.82T0.15. The difference in resulting model 
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Figure 4. From top to bottom are normalized STIS spectra of 
GJ86B obtained with the G430L, G430M, and G750M gratings. 
The only detected spectral features are the C2 Swan bands. Upper 
limits on calcium and hydrogen were derived from the Call K and 
Ha regions (Tabic [3j - 



parameters is significant and most likely caused by the two 
shortest wavelength fluxes. Indeed, the F225W and F275W 
bandpasses span a region influenced by at least one strong 
C 1 line which models tend to over-predict (in the core or 
wings ) . An example of this was observed bv lProvencal et all 
( 2002) in the STIS ultraviolet spectrum of Procyon B, where 
the strength of the 2478 Aline is much weaker (by nearly a 
factor of 30 in C/He) than predicted by models. The model 
shortcomings are understood to be a consequence of the 
high density environments in cool helium-rich white dwarf 
photospheres; the classical Van der Waal broadening treat- 
ment within the impact approximation, traditionally used 
in DQ model atmosp heres, is not appropriate for these lines 
i|Koester et al.ll 19821 ). 

If one ignores the two ultraviolet fluxes affected by these 
strong C 1 lines, and uses the remaining six WFC3 photo- 
metric data points, then T off = 8210 ± 180 K, log (C/He) = 
—4.76 ± 0.15 is obtained, in agreement with the param- 
eters from fitting the spectrum alone. For all stellar and 
evolutionary parameters, the weighted average of the pho- 
tometric (ignoring F225W and F275W) and spectroscopic 
effective temperatures, and the non-weighted average of the 
carbon abundances are thus adopted and listed in Table 
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A (A) A (A) 

Figure 5. Comparison between the observed and modeled stellar fluxes. In the left-hand panel, the model fluxes agree well with the 
WFC3 photometry in all bandpasses except F225W and F275W (see §3.1). In the right-hand panel the model spectrum reproduces both 
the optical continuum and C2 bands in the STIS G430L data, but indicates a strong C I 2478 A feature that under-predicts the stellar 
flux in the two ultraviolet filters. All spectral features are due to atomic or molecular carbon. 



(T cff = 8180 ± 120 K and log (C/He) = -4.8 ± 0.2). Upper 
limit hydrogen and calcium abundances were also derived 
and are listed in Table [3J while Figure [5] shows the best 
fitting model compared to the photometry and G430L spec- 
trum. All the derived stellar parameters are listed in Table 

El 



3.2 Current Binary Orbit 

The epoch 2012.2 WFC3 imaging provides a new point along 
the binary orbit of GJ 86, and is shown in F igure [ 6 ] together 
with previous astro metric measurements (|Lagrange et al.l 
l200fj ; lEls et~ai]|200ll ). Unfortunately, one quarter of the full 
orbit has not yet been observed, and hence the orbital pa- 
rameters can only be loosely constrained. For all orbital cal- 
culations, Ma = 0.8 Mq and Mb = 0.6 Mq are adopted 
for simplicity; the actual component masses are not better 
constrained at present. 

Using the new astrometric datum, a statistical anal- 
ysis of orbits that fit the observations was performed. A 
Monte Carlo simulation of 10 000 orbital solutions was cre- 
ated with periods from 20 to 500 yr, eccentricities from 0.00 
to 0.99, and with time of periastron passage from 1800 to 
2200. From these solutions, the 100 orbits with the best 
X 2 values, and with total masses within the limited range 
Mtot = 1.4 ± 0.1 Mq were selected. These orbits are plotted 
in Figure [5] and the resulting range in orbital parameters 
is listed in Table [4] As can be seen, only the system incli- 
nation is significantly constrained at present. The solutions 
with low eccentricities tend to have relatively short orbital 
periods, whereas the longer period solutions require high ec- 
centricities. Notably, the orbital solution with the shortest 
period yields a semimajor axis a = 27.8 AU with eccentric- 
ity e = 0.10, and is significantly larger than the representa- 
tive so lution with a = 18.4 AU adopted bv lLagrange et all 
<|200rj ). 



4 BINARY & PLANETARY EVOLUTION 
4.1 Total System Age 

The white dwarf cooling age places a fi rm lower limit 
of 1. 3 Gyr on the lifetime of the system (jFontaine et al.l 
2001). For the total age, one must add the hydrogen- 
burning lifetime of the progenitor, which should be in 
the range 1.2 — 1.6 Gyr for stars of mass 1.8 — 2.0 Mq. 
Thus, the evolved companion suggests a total system 
age of 2.65 ± 0.25 Gyr. Age estimates fo r the primary 
include 2.0 - 2.9 Gyr (ISaffe et all 120051 ) and 2.4 Gyr 
(|Mamaiek fc Hillenbrand! 120081 ) based on activity-age rela- 
tions, and these broadly agree wit h that derived for the 
secondary. iMamaiek fc Hillenbrand! (|2008l ) prefers an age 
of 3.7 Gyr based on activity-rotation-age metrics, but this 
value is significantly older than expected for the total life- 
time of the white dwarf and its progenitor. 



4.2 The Main-Sequence Progenitor System 

To constrain formation scenarios for the giant planet or- 
biting GJ 86 A, the binary orbit and component masses 
when both stars were on the main sequence is needed. Us- 
ing initial-to-final mass relations with some empirical con- 
straints at the low-mass end of progenito r and remnant 
masses (|Williams et al.l 120091 ; iKalirai et ai1l2008l ). GJ86B 
likely descended from a main-sequence, A-type star with a 
mass near 1.9 Mq. With both stars on the main sequence, 
the semimaj or axis of the binary was smaller by a factor 
(|jeanslll924f ) 



(1) 



oo _ M B + M A 
a ~ M Bo + M A 

This ratio is not very sensitive to the range of al- 
lowed masses, and yields ao = 0.52a for Ma = 
0.80M Q , M b = O.59M , M Bo = 1.86M Q . A canoni- 
cal mass value has been assumed for the primary, but 
is consist ent with isochrones a n d va ri ous literature es- 
timates (Ivan Belle fc von Braunl 120091: IValenti fc Fischerl 



|2005| ; iRibas et all 120031 ; iFlvnn fc Morelll Il997l ). From the 
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Table 4. Orbital Constraints from Monte Carlo Simulations 



Figure 6. The blue data points are the astrometric offsets de- 
termined he re using WFC3 ( leftm ost p oint), tog ether with those 
reported by lLagrange et al,l j2006T) and lEls et al, I 1 120011) . An ini- 
tial sample of 10 000 orbital solutions were generated through 
a Monte Carlo search, and overplotted in black are the 100 or- 
bits with the best \ 2 values, within the restricted mass range 
Mtot = 1.4±O.1M . 



current binary orbit simulations, a lower limit on the ini- 
tial semimajor axis is obtained ao 14.5 AU. Even at this 
lower limit, the 0.11 AU planetary orbit of G J 86Ab is stable 
for binary eccentricit ies e < 0.9 l|Holman fc Wiegertlll999l ; 
iRabl fc Dvora^ll988D . 



4.3 Planetesimals in GJ 86 

Although the class of detached binary stars with a white 
dwarf and a main-sequence component of spe ctral type K or 
earlie r are referred to as Sirius-type binaries (|Holberg et al.l 
2008), GJ 86 is also notably Procyon-like. Table[5]lists the at- 
mospheric and stellar parameters for these two white dwarf 
companion systems, demonstrating their broad similarities. 

Because GJ 86 is known to host a planetary system, 
with at least one giant planet orbiting the primary, it may 
also support a planetesimal population. These surviving mi- 
nor bodies could be gravitationally perturbed into close en- 
counters with the white dwarf, where they would be tidally 
shredded and eventua lly accreted , polluting the atmosphere 
with heavy element s J Jural 20031). as observed at nearly 30 



2012l ; lFarihi et alj|2012t 



single white dwa r fs (iGansicke et aL 

iKilic et all |2012| ; iDebes et al.l 120121 ). In a binary such as 
GJ 86, stable planetesimal belts can persist as circumstel- 
lar rings at either star, or a circumbinary disk. 

One region that remains stable over both the current 
and former binary configurations and e < 0.1 is within 
roughly 3AU of the K star primary. Stable circumstellar 
regions for the white dwarf are complicated by competing 
processes: while the region within 2 AU will be depleted of 
small bodie s by direct engulfment a nd tides during the gi- 
ant phases |Mustill fc Vil lavcr 201^), surviving objects will 



Parameter 


Range 


P 


120 - 481 yr 


To 


1933 - 2067 


e 


0.00 - 0.61 


a 


2'.'57 - &.'47 




27.8 - 69.8 AU 


i 


114? 7 - 122? 6 


n 


63? 7 - 76? 1 


Ul 


0° - 358° 



Note. The upper range of periods and semimajor axes result from 
the simulation input P m ax = 500 yr, but these orbits require high 
eccentricities e > 0.5 and are thus less likely. 



have their orbit expanded due to mass loss, by a factor 
of Mb /Mb ~ 3.2. However, the critical radius for stabil- 
ity also increases due to mass loss, but by a factor smaller 
than (M B + M A )/(M Bo + M A ) w 1.9, and therefore a nar- 
row range of stable, planetesimal orbits are possible for the 
white dwarf if a > 30 AU currently. It is worth noting that 
these calculations are purely static, and thus ignore the dy- 
namic effects of mass loss which may have significant im pact 
on stable orbits l|Veras fc Touul2012l ; rVeras et al.ll201ll ). Cir- 
cumbinary orbits beyond 80 AU are also stable in the current 
and former binary configurations for a « 30 AU and e < 0.1. 

Thus, one of the science goals of the HST observing 
program was to search for evidence of planetesimals in 
this system, via atmospheric metal pollution in the white 
dwarf. As discussed in §2.2, calcium absorption was not 
detected in the STIS data, and a strict upper limit of 
log (Ca/He) < —11.8 was derived for GJ86B. This corre- 
sponds to Men < 3.2 x 10 17 g in the stellar convection zone, 
or an upper limit to the total accreted heavy element mass 
of M z < 2.0 x 10 19 g for material that is 1.6 % calcium by 
mass as in the bulk Earth l|Allegre et al.lll995l ). Because the 
Call K line is the most sensitive tracer of metal pollution 
at optical wavelengths for whi t e dwarfs of this effective tem- 
perature (jKoester et al.l [20051 ; IZuckerman et alj |2003). this 
large cometary-sized mass is the strongest possible limit to 
the observable planetesimal population. 



4.4 A Failed Planetary System at Procyon? 

Interestingly, and also using STIS, iProvencal et al.l (|2002T ) 
found that Procyon B exhibits a strong Mg II resonance line 
at 2800 A, a few Fe lines in the near-ultraviolet, and Call 
H and K absorption in the optical. While Procyon is not 
known to host any planets or planetesimal debris (via pre- 
cision radial velocity, transits, astrometry, direct imaging, 
or infrared excess), the white dwarf is externally polluted 
with at least three heavy element s that could have origi- 
nated in a minor planetary body (|Farihi et al. 2010). Us- 
ing the published magnesium abundance, Procyon B has ac- 
creted at least Myig = 5.1 x 10 18 g within the past few 10 7 yr, 
where 2.8 Myr is a single sinking ti mescale for th is element 
at the stellar effective temperature (|Ko cstcr 2009]). For ma- 
terial with a chondritic or Earth-like elemental abundance, 
this is equivalent to a total accreted heavy element mass of 
at least 3 - 6 x 10 19 g (|Lodders fc Feglevll 19981 ; lAllegre et all 
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Table 5. Comparison of GJ 86B & ProcyonB 



PeLrfiii'H 'tor 


GT86B 


Procyon B a 


SpT 


DQ6.2 


DQZ6.5 


Tefi (K) 


8180 


7740 


Mass (Mq) 


0.59 


0.60 


log(H/He) 


< -4.3 


< -3.7 


log(C/He) 


-4.8 


-5.5 


log (Mg/He) 6 




-10.4 


log (Ca/He) 


< -11.8 


-11.8 


Primary SpT 


KOV 


F5IV 


a (AU) 


> 28 


15.2 


e 




0.41 



a Atmospheric parameters from [Provencal et al.l l|2002h . 
b The STIS spectra of GJ 86B do not cover the Mgn 2800 A res- 
onance lines. 

1 19951 ). and equivalent to a 25 — 35 km diameter object for 
typical asteroid densities. 

Other possible sources for the atmospheric metals are 
the interstellar medium , or the stellar wind from Procyon A 
|Provencal et alj|2002j ). In both of these cases the accreted 
matter would have solar abundances (lo g (Mg/H) = —4.45, 
log (Mg/He) = -3.35; lLoddersI 120031 ) and imply time- 
averaged accretion rates of Mh = 8.0 x 10 8 g s" 1 and MHe = 
2.6 x 10 8 gs -1 . If the heavy elements have been continually 
captured from Procyon A, then hydrogen and helium should 
have accumulated i n the white dwarf at mosphere over its 
1.4 Gyr cooling age |Fontaine et all [200 if ). During this pe- 
riod 3.6 x 10 25 g of hydrogen and 1.1 x 10 25 g of helium will 
have been added to the 2.1 x 10 28 g helium convection zone, 
yielding a current abundance log(H/He) = —2.2. However, 
the upper limit derived from the lack of Ha absorption is 
log (H/He) < —3.7 and inconsistent with stellar wind accre- 
tion. 

If one ignores the issue of hydrogen accumulation, the 
mass loss rate from Procyon A that is necessary to supply 
the magnesium in the white dwarf at the observed time- 
averaged rate is 1.8 x 10~ 9 Mq yr" 1 for a Bondi-Hoyle accre- 
tion flow l |Debej2006h flow and 9.1 x 10" 8 Mq yr -1 for grav- 
itational (Eddington) accretion. These mass-loss rates are 
between 5 and 7 orders of ma gnitude larger than that mea- 
sured for the Sun and aCen l|Wood et al.ll200ll ; IWithbroe] 
1989). Therefore, wind capture can be ruled out with con- 
fidence. Similar arguments can be used to discount the ac- 
cretion of interstellar matter. First, hydrogen accumulates 
whereas metals sink, and thus the hydrogen-deficient atmo- 
sphere of ProcyonB is at odds with the accretion of solar 
abundance material for more than 2% of its cooling lifetime. 
Second, molecular clo ud densities are req uired to account for 
the observed metals (|Farihi et alj 12010 ). but Procyon has 
been moving within the l ocal, r ~ 100 pc, ISM-poor Bubble 
l|Redfield fc Linsky|[200ct) for at least 2.6 Myr. 

The remaining possibility is the accretion of one or more 
planetesimals. Accounting for both the current and former 
binary configurations, stable orbital regions at Procyon can 
be either circu mbinary beyond 55 . 4 AU or circumprimary 
within 2.3 AU (|Holman fc Wiegertl I1999T ) . although bodies 
in the latter region would be dynamically unlikely to be per- 
turbed towards the surface of the white dwarf. There is no 



circumstellar region at ProcyonB that remains stable over 
the lifetime of the binary. Thus, a circumbinary population 
of planetesimals is the most likely source of the metals ob- 
served in ProcyonB. Observations with Herschel PACS at 
160 fim have not detected a far-infrared excess at Procyon, 
and limit the fractional dust luminosity to L/L* < 5 x 10~ 7 
over the 30 — 300 AU region (G. Kennedy 2012, private com- 
munication). This limit is only a few times gre ater than the 
fract i onal dust lu minosity of the Kuiper belt (| Booth et al.1 
l200gt ISternlll996ft . and for similar particles would result in 
an upper dust mass limit of roughly an Earth mass for Pro- 
cyon. 

On the main sequence, Procyon should have had a semi- 
major axis near 9 AU. Even for e = 0.1 this implies that disk 
material available for planet building would be confined to 
within 2.0 AU of Procyon A and 2.4 AU of ProcyonB, and 
thus completely within their respective snow lines. Thus, gi- 
ant planet formation at either component of the Procyon 
system is likely precluded, and likely also for the a > 55 AU 
stable circumbinary environment. Only the inner regions of 
Procyon A appe ar capable of forming small, solid planets 
(as in q Cen AB; iDumuscme et al.l l2012). that persist to the 
present day, assuming the binary eccentricity was not larger 
when both stars were on the main sequence. Thus, the Pro- 
cyon system may represent a case of failed or truncated 
planet formation, where large planetesimals were formed but 
further growth was prohibited. 



5 SUMMARY & OUTLOOK 

High-contrast optical imaging and low-resolution spec- 
troscopy with HST have unambiguously characterized the 
stellar companion to the planet-host star GJ 86A. The sec- 
ondary is a relatively cool white dwarf with a helium- 
dominated atmosphere and molecular absorption bands due 
to trace carbon, which is almost certainly dredged from 
the core. Neither Ha nor Call K are detected in deeper, 
medium-resolution spectroscopy, the latter placing modest 
limits on the presence of scattered planetesimal material. 

The binary separation of GJ 86 has continued to in- 
crease since last observed in 2005, and the actual orbit de- 
viates from prior estimates. While one-quarter orbit has not 
yet been observed, Monte Carlo simulations constrain the 
semimajor axis to a ^ 27.8 AU, and the system inclination 
within 114?7 ^ i 122?6. For planetary-binary coplanar 
orbits, these results imply the mass of GJ 86Ab lies between 
4.4 and 4.7M Jup . 

Future imaging is necessary to observe apastron and dy- 
namically constrain the component masses. This can be done 
with adaptive optics in the near-infrared, but superior data 
will come from follow up imaging and spectroscopy with 
HST in the near-ultraviolet. In this relatively low-contrast 
wavelength regime, there will be narrow absorption lines 
from C I, and possibly Mg n which is a strong absorption fea- 
ture in ProcyonB. Such spectral lines can provide a baseline 
radial velocity and gravitational redshift for the secondary; 
these would more efficiently characterize the current and 
former binary configurations than near-infrared data. 



© 0000 RAS, MNRAS 000, 000-000 



Constraints on Planet Hosting Binary GJ 86 9 



ACKNOWLEDGMENTS 

The authors thank an anonymous referee for a timely and 
thorough report. J. Farihi thanks J. Provencal for helpful 
feedback on the available spectra of Procyon B, as well as G. 
Kennedy and M. Wyatt for a useful summary of the infrared 
excess and dust limits in this system. This work is based on 
observations made with the Hubble Space Telescope which is 
operated by the Association of Universities for Research in 
Astronomy under NASA contract NAS 5-26555. These ob- 
servations are associated with program 12548. Support for 
Program number 102548 was provided by NASA through 
grant HST-GO-12548 from the Space Telescope Science In- 
stitute. J. Farihi gratefully acknowledges the support of the 
STFC via an Ernest Rutherford Fellowship and as a PDRA. 
N. Haghighipour acknowledges support from NASA grants 
EXOB NNX09AN05G and HST-GO-12548.06-A. J. Holberg 
acknowledges support from NSF Grant AST-1008845. M. 
Barstow and M. Burleigh acknowledge support from the 
STFC. 



REFERENCES 

Akeson R. L., Koerner D. W., Jensen E. L. N. 1998, ApJ, 
505, 358 

Allegre C. J., Poirier J. P., Humler E., Hofmann A.W. 1995, 

Earth Planetary Sci. Letters, 4, 515 
Artymowicz P., Lubow S. H. 1994, ApJ, 421, 651 
Booth M., Wyatt M. C, Morbidelli A., Moro-Martm A., 

Levison H. 2009, MNRAS, 399, 385 
Correia A. C. M., et al. 2008, A&A, 479, 271 
Debes J. H. 2006, ApJ, 652, 636 

Debes J. H., Kilic M., Faedi, F., Shkolnik E. L., Lopez- 
Morales M., Weinberger A. J., Slesnick C, West R. G. 
2012, ApJ, 754, 59 
Desidera S., Barbieri M. 2007, A&A, 462, 345 
Dressel L., et al. 2010, Wide Field Camera 3 Instrument 

Handbook, Version 3.0, (Baltimore: STScI) 
Dufour P., Bergeron P., Fontaine G. 2005, ApJ, 627, 404 
Dumusque X., et al. 2012, Nature, 491, 207 
Els S. G., Sterzik M. F., Marchis, F., Pantin E., Endl M., 

Kiirster M. 2001, A&A, 370, LI 
Farihi J., Barstow M. A., Redfield S., Dufour P., Hambly 

N. C. 2010, MNRAS, 404, 2123 
Farihi J., Gansicke B. T., Steele P. R., Girven J., Burleigh 

M. R., Breedt E., Koester D. 2012, MNRAS, 421, 1635 
Flynn C, Morell O. 1997, MNRAS, 286, 617 
Fontaine G., Brassard P., Bergeron P. 2001, PASP, 113, 
409 

Gansicke B. T., Koester D., Farihi J., Girven J., Parsons 

S. G., Breedt E. 2012, MNRAS, 424, 333 
Ghezzi L., Cunha K., Smith V. V., de Araujo F. X., Schuler 

S. C, de la Reza R. 2010, ApJ, 720, 1290 
Haghighipour N. 2006, ApJ, 644, 543 

Hatzes A. P., Cochran W. D., Endl M., McArthur B., Paul- 
son D. B., Walker G. A. H., Campbell B., Yang S. 2003, 
ApJ, 599, 1383 

Hayashi C. 1981, Prog. Theor. Phys. Suppl., 70, 35 

Holberg J. B., Sion E. M., Oswalt T., McCook G. P., Foran 
S., Subasavage J. P. 2008, AJ, 135, 1225 

Holman M. J., Wiegert P. A. 1999, AJ, 117, 621 



Jeans J. H. 1924, MNRAS, 85, 2 
Jura M. 2003, ApJ, 584, L91 

Kalirai J. S., Hansen B. M. S., Kelson D. D., Reitzel D. B., 
Rich R. M., Richer H. B. 2008, ApJ, 676, 594 

Kilic M. Patterson A. J., Barber S., Leggett S. K., Dufour 
P. 2012, MNRAS, 419, L59 

Koester D. 2009, A&A, 498, 517 

Koester D., Rollenhagen K., Napiwotzki R., Voss B., 
Christlieb N., Homeier D., Reimers D. 2005a, A&A, 432, 
1025 

Koester D., Weidemann V., Zeidler E. M. 1982, A&A, 116, 
147 

Lagrange A. M., Beust H., Udry S., Chauvin G., Mayor M. 

2006, A&A, 459, 955 
Lodders K. 2003, ApJ, 591, 1220 

Lodders K., Fegley B. 1998, The Planetary Scientist's Com- 
panion (New York: Oxford University Press) 
Mamajek E. E., Hillenbrand L. A. 2008, ApJ, 687, 1264 
Mathieu R. D. 1994, ARA&A, 32, 465 
McCook G. P., Sion E. M. 1999, ApJS, 121, 1 
Mugrauer M., Neuh&user R. 2005, MNRAS, 361, L15 
Mustill A. J., Villaver E. 2012, ApJ, 761, 121 
Muterspaugh M. W., et al. 2010, AJ, 140, 1657 
Nelson A. F. 2000, ApJ, 537, L65 

Pelletier C, Fontaine G., Wesemael, F., Michaud G., Weg- 

ner G. 1986, ApJ, 307, 242 
Perryman M. A. C, et al. 1997, A&A, 323, L49 
Prato L., Weinberger A. J. 2010, Planets in Binary Star 

Systems, ed. N. Haghighipour, Springer, New York, 1 
Provencal J. L., Shipman H. L., Koester D., Wesemael, F., 

Bergeron P. 2002, ApJ, 568, 324 
Queloz D., et al. 2000, A&A, 354, 99 
Rabl G., Dvorak, R. 1988, A&A, 191, 385 
Raghavan D., Henry T. J., Mason B. D., Subasavage J. P., 

Jao W. C, Beaulieu T. D., Hambly N. C. 2006, ApJ, 646, 

523 

Redfield S., Linsky J. L. 2008, ApJ, 673, 283 
Ribas I., Solano E., Masana E., Gimenez A. 2003, A&A, 
411, L501 

Rodriguez L. F., et al. 1998, Nature, 395, 355 

Saffe C, Gomez M., Chavero C. 2005, A&A, 443, 609 

Stern S. A. 1996, A&A, 310, 999 

Valenti J. A., Fischer D. A. 2005, ApJS, 159, 141 

van Belle G. T., von Braun K. 2009, ApJ, 694, 1085 

van Leeuwen F. 2007, A&A, 474, 653 

Veras D., Tout C. A. 2012, MNRAS, 422, 1648 

Veras D., Wyatt M. C, Mustill A. J., Bonsor A., Eldridge 

J. J. 2011, MNRAS, 417, 2104 
Weidenschilling S. J. 1977, ApSS 51, 153 
Williams K. A., Bolte M., Koester D. 2009, ApJ, 693, 355 
Withbroe G. L. 1989, ApJ, 337, L49 

Wood B. E., Linsky J. L., Miiller H. R., Zank G. P. 2001, 
ApJ, 547, L49 

Zucker S., Mazeh T., Santos N. C, Udry S., Mayor M. 

2004, A&A, 426, 695 
Zuckerman B. Koester D., Reid I. N., Hiinsch M. 2003, 

ApJ, 596, 477 



© 0000 RAS, MNRAS 000, 000-000 



